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Actuator rate saturation is an important factor adversely affecting the stability and performance of aircraft 
flight control systems. It has been identified as a catalyst in pilot-induced oscillations, some of which have been 
catastrophic. A simple design technique is described that utilizes software rate limiters to improve the performance 
of control systems operating in the presence of actuator rate saturation. As described, the technique requires control 
effectors to be ganged such that any effector is driven by only a single compensated error signal. Using an analysis 
of the steady-state behavior of the system, requirements are placed upon the type of the loop transmissions and 
compensators in the proposed technique. Application of the technique to the design of a multi-input/multi-output, 
lateral-directional control system for a simple model of a high-performance fighter is demonstrated as are the 
stability and performance improvements that can accrue with the technique. 


I. Introduction 

T HE performance requirements of modem, high-performance, 
fly-by-wire aircraft have made it imperative that the charac- 
teristics of the actuation devices be included in any control system 
design procedure. For example, in the control of supermaneuverable 
aircraft, the reality of control actuator limitations can determine the 
overall control system design philosophy. 1 Although the linear dy- 
namics of the actuator are often modeled in such designs, the non- 
linear behavior, e.g., saturation, is not often considered in explicit 
fashion. Typically, extensive a posteriori simulation is employed 
to determine if and when actuator saturation is likely to occur and 
then if such saturation can affect flight safety. Even time-consuming 
simulation efforts are often unable to uncover situations that can af- 
fect flight safety. For example, probably no aircraft in the history of 
aviation has undergone more simulator evaluation than the Shuttle 
Orbiter. Yet a serious pilot-induced oscillation (PIO) occurred in 
early free-flight testing, 2 which led to actuator rate saturation. Ac- 
tuator rate saturation has also been implicated in the YF-22 crash 
early in 1992 (Ref. 3) and possibly in the JAS 39 crash of 1993 (Ref. 
4), both of which also experienced serious PIOs. Indeed, the pivotal 
role that actuator rate saturation may play in PIOs is coming under 
increased scrutiny. 5 

II. Background 

Feedback design techniques that explicitly consider the possi- 
bility of control saturation are receiving increasing attention in the 
literature. Optimal control strategies for regulators employing in- 
equality constraints on control and state variables are well known. 6 
Techniques that attempt to minimize the adverse effects of saturation 
are also well established, e.g., antiwindup controllers. 7 Feedback 
systems designed to avoid saturation have been discussed. 8 Nonlin- 
ear feedback controllers designed with the goal of improved perfor- 
mance with state and control saturation have also been described. 9 
Linear feedback schemes to improve the performance of systems 
with saturation nonlinearities have been formulated. 10 Finally, an 
alternate control scheme has been recently proposed to reduce the 
adverse effects of time delays associated with flight control system 
actuator rate limiting, particularly that which occurs in pilot-induced 
oscillations. 5,11 

As opposed to the methods just mentioned, the research to be de- 
scribed employs a technique for addressing actuator rate saturation 


Received Feb. 12, 1996; revision received Aug. 6, 1996; accepted for 
publication Oct. 18, 1996. Copyright © 1996 by the American Institute of 
Aeronautics and Astronautics, Inc. All rights reserved. 

‘Professor, Department of Mechanical and Aeronautical Engineering. As- 
sociate Fellow AJLAA. 

t Assistant Professor, Department of Mechanical and Aeronautical Engi- 
neering. Member AIAA. 


that is very simple in theory and implementation. Closed-loop sta- 
bility and performance are improved through the introduction of 
software rate limits, i.e., nonlinear systems that ensure that no rate 
commands are sent to the actuators that are beyond the capabili- 
ties of these devices. Software rate limits are certainly not novel; 
however, it is believed the technique to be described uses them in a 
unique manner. As will be seen, although the proposed design tech- 
nique is very simple, it does require a particular feedback topology 
for its implementation. 

Section HI introduces the technique using a single-input/single- 
output (SISO) formulation and example. Section IV extends the 
technique to multi-input/multi-output (MIMO) designs and demon- 
strates the aforementioned limitations on feedback topology. This 
section also applies the technique to a lateral-directional flight con- 
trol problem using a simple model of a high-performance fighter 
aircraft. Finally, a brief discussion of the design technique is pre- 
sented in Sec. V, and conclusions are drawn in Sec. VI. 

HI. Software Rate Limits: SISO Systems 
Software Limiter 

Consider Fig. la, which shows a very simple SISO control sys- 
tem. Although the linear actuator dynamics are ignored here, it will 
be assumed that actuator rate limiting does occur. The rate limited 
actuator A RL is modeled as shown in Fig. lb. Figure 2 shows the 
proposed software limiter S RL as part of a SISO feedback system. 
S RL consists of a derivative element, followed by an amplitude lim- 
iter, followed by an integration element As is evident from their 
structures, the input-output behavior of S R l and A RL (or an actual 
rate-limiting actuator) will be different. The former device employs 
no internal feedback, whereas the latter does. Consider the behav- 
ior of each element in isolation as shown in Fig. 3. The input is a 
ramp followed by a constant value. Although both elements result in 
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"software" limiter, S RL 



Fig. 2 System of Fig. la with a software rate limiter. 



the assumption of closed-loop stability [E(s) with no right half- 
plane poles]. Thus, invoking the final value theorem 

e ss (t ) = lim e(t ) = lim sE (3) 

t~* oo S-+0 


But, from Eq. (2), 


sE = [sR/(\+Ly\-D n P (4) 

Now denote the type of a transformed variable or transfer function A 
as q A . Type refers to the exponent on any free s in the denominator 
of die transformed variable or transfer function. 12 Then the types of 
the first and second terms on the right-hand side of Eq. (4) are 


Fig. 3 Rate limiting behavior: a) input/output behavior of software 
limiter of Fig. 2 and b) input/output behavior of rate limiter of Figs. 1 
and 2. 

rate limited outputs, these outputs are quite different. The software 
limiter in Fig. 3a merely prevents rates beyond a certain magni- 
tude from being sent to the actuator. It comes out of saturation as 
soon as |x(/)| < M. The output of the rate limit element of Fig. lb 
(and a rate limiting actuator) stays in saturation until the displace- 
ment error signal E v becomes appropriately small. It is this latter 
behavior that can seriously compromise performance and stability 
when the actuator is part of a closed-loop system. As will be seen, 
the performance of a system with the software limiter and a rate- 
limiting actuator can be far superior to a system without the soft- 
ware limiter. However, the behavior of the software limiter leads to 
requirements upon the system loop transmission if steady-state per- 
formance is to be maintained. This will be demonstrated in the next 
section. 

Quasilinear System and Steady-State Error 

Consider Fig. 4, a quasilinear representation of Fig. 2, in which 
the limiter has been replaced by a remnant signal d(t). The remnant 
is defined simply as that signal that must be added to x(t) in Fig. 4 
so that y(t) in Fig. 4 is identical to y(t) in Fig. 2 when limiting 
occurs. In what follows, the dynamics of the rate limiter A RL can 
be neglected, because saturation of A RL will be eliminated by Sri 
and the linear dynamics of A RL are negligible (k 2> 1 in Fig. lb). 
Block diagram algebra can then be used to show the following: 

Y = D + sGR-(sGPY/s) = D + sGR-LY (la) 
(1 + L)y = D+sGR (lb) 

D/{ 1 + L) = D n = Y — [sGR/(l + L)) = Y~Y l (lc) 


q Uttena = q R - q L-l (5a) 

q 2x<Ucrm = q D n+q P (5b) 

where it has been assumed in obtaining Eq. (5a) that q L > 0, a 
criterion invariably met in any satisfactory feedback design. Now, 
unfortunately, D n cannot be treated like an independent disturbance 
as evident from Eq. (lc), i.e., it will be dependent on the characteris- 
tics of the system into which it is injected, e.g., the loop transmission 
L. If q L > q p , then a necessary and sufficient condition for the sta- 
ble linear system of Fig. 2 to exhibit asymptotic stability in response 
to d(t) is q Dn < —q L . A conservative expression for q Dn could then 
be offered as 

q Dn = ~q L (6) 

Using Eq. (6), the right-hand side of Eqs. (5a) and (5b) will be 
negative and e ss (t) -*• 0 if 

From Eq. (5a): 

q L > q R (7a) 

From Eq. (5b): 

q L > q p + 1 or q G > 1 (7b) 

Although dependent on the assumption of stability, inequalities 
(7) provide design guidance and offer necessary conditions for 
asymptotic stability. The requirements implied by inequalities (7) 
are necessitated by the aforementioned behavior of the software 
limiter, i.e., y(t) not approaching x(t) after a period of limiting. In 
general, the larger the*? 1 , the better the performance of the system 
with the software limiter, provided that adequate stability margins 
can be maintained in the linear design. 


where Y L is the signal that would exist in either Fig. 2 or Fig. 4 
when no limiting occurs. Note from Eq. (1) that D n is merely the 
difference between Y (with limiting) and Y L . 

Now, from Fig. 4 

E = [/?/(l + £)] — (D n P/s) (2) 

A fundamental assumption in this analysis is that the closed-loop 
system of Fig. 2 is stable, even in saturation. This will allow the final 
value theorem to be applied to Eq. (2). Note that open-loop stability 
is not required, because if P(s) has right half-plane poles, they will 
be canceled by identical poles of D n (s), a condition guaranteed by 


SISO Example 

Consider the system of Fig. 1 . Let 


The compensator G has been designed so that the loop transmission 
has a crossover frequency of 2 rad/s, with desirable loop transmis- 
sion characteristics. The amount of rate limiting will be adjusted to 
produce significant saturation in this example. Assuming that 

q* < 1 (9) 
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Fig. 4 Quasiiinear representation of Fig. 2. 




Fig. 5 Response of system of Fig. la to pulsive input, with and without 
rate saturation. 


Fig. 6 Response of system of Fig. 2 to pulsive input, with and without 
rate saturation. 


then inequalities (7) require that 

q L > q p + 1 = 2 (10) 

Inequality (10) requires that the type of the loop transmission L = 
G P be increased from the type 1 implied by Eq. (8) when the system 
of Fig. 2 is implemented. Thus, one can define a new compensator as 


s 

This modification of the original compensator G will cause some 
change in the performance of the system of Fig. 2 as compared with 
that of the system of Fig. la when no saturation occurs; i.e., there 
will be a modest reduction in stability margins. Finally, the s appear- 
ing before the limiter in Fig. 2 will be included in the compensator 
in actual implementation. Although the resulting compensator sG' 
is now only proper, as opposed to being strictly proper, the integra- 
tor following the limiter will restore the latter property. This can 
be important from the standpoint of sensor noise propagation to the 
plant actuator. 

Figure 5 shows the response of the system of Fig. la to a pulsive 
input of unity magnitude and 5-s duration with and without rate satu- 
ration. For the case with saturation, the rate limit was 0.2/s. Figure 6 
shows the response of the system of Fig. 2 [with compensator G'{s) 
and the software limiter] to the same input, with and without the 
rate limit of 0.2/s. The performance improvement with saturation is 
obvious. An examination of the actuator rate time histories in this 
example indicated that the system of Fig. 2 was in saturation far 
less than the system of Fig. la. A comparison of the dashed curves 
in Figs. 5 and 6 show that the performance of the system of Fig. 2 
is somewhat poorer than that of Fig. la when no saturation occurs. 
This is attributable to a slight reduction in stability margins when 
q L was increased from 1 to 2 [using G'(s) as opposed to G(s)] as 
dictated by inequality (10). However, this is a small price to pay for 
the performance improvement that occurs with saturation. 

Figure 7 compares the signals u(t ) for the systems of Figs. 1 a and 
2 when a sinusoidal input 

( 12 ) 



Fig. 7 Comparison of signals u(t) in systems of Figs, la and 2 with 
sinusoidal input and rate saturation. 

was applied. Here the rate limit was set to 2.0/s. The characteristics 
of u(t) indicate that both systems were in nearly constant rate satu- 
ration of ±2.0/s throughout. Because of this significant limiting, the 
tracking behavior of both the systems of Figs, la and 2 was poor. 
What is important, however, is the fact that the response u(t) of 
Fig. la exhibits an effective incremental time delay of 0.55 s com- 
pared with the same system without rate limits. In comparison, the 
response u(t) of Fig. 2 exhibits an effective incremental delay of 
only 0.12 s compared with the same system without rate limits. The 
effective incremental delays were calculated between correspond- 
ing instants when u(t ) changes sign. If the systems of Figs, la and 2 
were representing the stability and command augmentation system 
of an aircraft, with r(t) being the pilot’s input, the much larger delay 
of the system of Fig. la would probably lead to this system being 
more prone to PIOs than the system of Fig. 2 (Ref. 5). 

IV. Software Rate Limits: MEMO Systems 
Software Limiters 

Any application of the technique introduced in Sec. m to a re- 
alistic aircraft flight control problem must allow consideration of 


r(t) = sin(3r) 
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MIMO systems. This section is directed to that end. For purposes 
of exposition, a 2 x 2 system will be employed, i.e., two inputs and 
two controlled outputs. Figure 8a shows the nonlinear system, and 
Fig. 8b shows its quasilinear representation. Note that n e control 
effectors are allowed (where n e can be any number). However, the 
effectors must be ganged; i.e., a group of effectors or actuators is 
driven by only one of the n r compensated error signals. Here n r — 2. 
In this case, the ganging is done by matrix K, each row of which 
contains only one nonzero entry. These elements are numerically 
equal to the rate limit of the ith actuator. Thus, in Figs. 8a and 8b, 
none of the S Cj will exceed the rate limits of the associated actuators. 
The K matrix can be succinctly expressed as 

K = R l ■ T (13) 

where R L is an n e x n e diagonal matrix with the actuator rate limits 
appearing along the diagonal and r is the n e x n r control distribution 
matrix. Each row of T contains only one nonzero element, which 
is unity. The justification for this choice of K will be presented in 
what follows. 

In Figs. 8a and 8b, the dynamics of the actuators (with limiting) 
are assumed to be included in the plant element Pij . Finally, one can 
define an effective plant matrix P e as 

P e = PK (14) 


limiters alone. Thus, the effects of actuator amplitude limiting can- 
not be considered in this analysis. 

The loop transmissions P, obtained by successively opening the 
loops in Fig. 8b at U ci are given as 


L\ = G\ 


Pe n + G 2 A 
1 + Perfil 


t ^ P e 22 + GlA 

i2 = G2 'TT^G7 


A = P P — P P 

U r <U r «22 r *12 r ‘2l 


By a procedure analogous to that used to obtain Eq. (4), one can 
show 

sRi sP en G 2 R 2 


l+^i (l+P^l + P^G,) 


All Pi _ Ai2-Pei 2 (l + A) 1 

G i (1 + Pi) (l + P en G%) 

__ s R 2 sP^GiRi 

~ I+P 2 (l+P 2 )(l + P en Gi) 


Quasilinear System and Steady-State Error ~ , n p n 4 - f 1 1 

Following a procedure analogous to that which yielded Eqs. (1), 2L_! /lL * 21 ^ — (17b) 

D ni and P>„ 2 in Fig. 8b can be shown to be G 2 (1 + P 2 ) (l + P en Gi) 


D ni = A/d+Pi) = Yi-Y lL D ni = D 2 /(l+P 2 ) = Y 2 -Y 2l 

(15) 

Each Y, l in Eqs. (15) can be defined as the signal Y-, that would exist 
if no rate limiting was occurring in the ith software limiter, alone. 
As Fig. 8a indicates, amplitude, as well as rate limiting, can occur 
in the control system actuators. However, the D ni in the quasilinear 
systems of Figs. 4 and 8b have been assumed to arise from software 


Each of the four terms on the right-hand sides of Eqs. (17a) and 
(17b) represents contributions from the four exogenous inputs Pi, 
R 2 , D „ j , and D„ 2 (remembering that D„ t are not really exogenous). 

Assuming that the closed-loop system of Fig. 8 a is stable, even 
with saturation, the final value theorem can be invoked as with the 
previous SISO system. As was done in Eqs. (5), one can determine 
the type of each of the terms on the right-hand side of Eqs. (17a) 
and (17b) as follows. 



Fig. 8 A 2 x 2 MEMO system: a) with software rate limiting and b) a quasilinear representation of panel a. 
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Equation (17a): 

q 1st term _ gRi _ qL\ _ j 
q2nd\erm _ qR 2 _j_ q p ‘n — q — ^^22 — 1 
^3rdterm _ qD„ x _j_ qL\ _ qG\ 

^ 4th term _ q D n 2 _j_ q p t\ 2 — q^l -j- q L 2 — q p t 22 — (JI G 2 

Equation (17b): 

qUtterm = qR 2 _ qL 2 __ j 

^ 2nd term __ ^/?i _j_ ^< 21 — q L2 — £^ e ii — 1 

^3nJterm _ ^£>« 2 _j_ ^G 2 

^ 4th term _ qD n[ _j_ qP t2l _ ^ L 2 q. _ qGi 


(18) 


(19) 


In obtaining Eqs. (18) and (19) it has been assumed that q L ‘ > 0 
and q p ‘u Gi > 0, criteria invariably met in the satisfactory design of 
the feedback system of Fig. 8a. The design task is to select q L ‘ and 
q Gi (i = 1, 2) so that 


qjlh term < _j (y = 1, 2, . . .) 

in Eqs. (18) and (19). As in Eqs. (6) and (9), allowing 

q D "i = -q L ‘ q R i < 1 

Eqs. (20) and (21) lead to Eqs. (18) and (19) requiring 
q Ll >1 q Ll > q P ‘ n — q Pt22 + 1 q Gi > 1 
q tl > q p 'tt — q Pe & — q° 2 + 1 


( 20 ) 

( 21 ) 


( 22 ) 


and 

q Li > 1 


q L 2 > q p “ 2 1 — + 1 q° 2 > 1 


> q Pe 2\ — + 1 


(23) 


The final value theorem and inequality (20) will ensure 


linear model - — (s) = 


\_k 

m 2 

+ 

h/ 
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Fig. 9 Model of a second-order amplitude and rate limited actuator. 


thrust vector input Sytv were used to control sideslip p, and the dif- 
ferential elevator input S DT , aileron input S A , and differential pitch 
thrust vector input S RT v were used to control roll rate p s . Obviously, 
other choices could have been made. In terms of Eqs. (13), 


"60 

0 

0 

0 

0 " 



"0 

1" 

0 

100 

0 

0 

0 



0 

1 

0 

0 

100 

0 

0 

deg/s 

r = 

i 

0 

0 

0 

0 

60 

0 



0 

1 

_ 0 

0 

0 

0 

60. 



_1 

0_ 


(25) 


e iss (t) = 0 i = l,2 (24) 

The q L ‘ are governed by the q G ‘ as evident in Eqs. (16a) and (16b). 
Of course, the methodology just outlined can be extended to 3 x 3, 
4x4 systems, etc. As in the case of the SISO design, open-loop 
stability is not required. 


k = r l • r = 


0 

60 " 

0 

100 

100 

0 

0 

60 

60 

0 


Ganged Controls 

The explanation for imposing ganged controls can now be given. 
First, grouping each software limiter with a compensator ensures 
that the type of each of the n r loop transmissions q Li (n r = number 
of controlled outputs) can be determined by a single compensator; 
e.g., in Eq. (16), q Ll can be governed by G lt etc. Second, ganging 
the effectors and defining the nonzero elements of the K matrix as 
in Eq. (13) allows n r software limiters to prevent rate saturation of 
n e actuators, with n e > n r . 


and 

8 = [ 5 or » 8 A , Sr, Srtv, 8 ytv] t 

A classical loop-shaping procedure was used to determine the 
compensators Gp and G p , shown in Table 1. The crossover fre- 
quencies of each loop were selected as 

co Cf> = 5.0 rad/s co Cp = 4.0 rad/s (26) 


MIMO Example 

The aircraft example to be examined involves the lateral- 
directional control of a model of a high-performance fighter air- 
craft. The simplified aircraft and actuator dynamics were taken from 
Ref. 13 and are repeated here in Table 1 . Rather than employing the 
simple rate limiter shown in Fig. lb, a more accurate representation 
of a rate and amplitude limited second-order actuator was employed. 
This actuator model is shown in Fig. 9 (Ref. 14). The controlled out- 
puts are sideslip P and roll rate p s (r) about the x stability axis. Note 
here that n r = 2 and n e = 5. Thus, in terms of Fig. 8b, Ri = p c Fp, 
Ci = P, R 2 = p c F p , C 2 — p s , G\ = Gp, and G 2 — G p . Selection 
of which of the five control effectors to associate with each of the 
compensators can be based upon the control axis authority of the ef- 
fectors themselves. Here, for example, the rudder input S R and yaw 


In this application, q P ^ = 0, and inequalities (22) and (23) merely 
require that 

qh > 1, qh > 1, q°i > 1, q <h > 1 (27) 

With the Gr and G p of Table 1, inequalities (27) would already be 
met with q^ = q L <> = 1. However, the types of G'^ and G' p were 
chosen so as to yield 


ih — q L P = 2 


(28) 


Thus, 


G'p^Gp 


(f + 0.5) 




(s + 0.4) 


(29) 
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Table 1 Vehicle dynamics and compensation 
Flight condition 

Altitude = 20,000 ft, M — 0.5, atrira = 5.4 deg 
Vehicle dynamics 

x=Ax+Bu X = [p,p,r] T 
u — [Sot, Sa, Sr, Srtv, <$trv] r 


" -0.1354 

0.09036 

-0.9949” 

A = -10.37 

-1.469 

0.5126 

2.181 

-0.01482 

-0.1277 

“ -0.01091 -0.005695 

0.01555 

0 0.00489 

B = 9.93 12.12 

0.9416 

0.3977 0.02817 

0.2757 -0.2797 

-0.7419 

-0.001175 -0.3861 

where 




State variables 

— sideslip angle, deg 
p = roll rate, deg/s 
r — yaw rate, deg/s 
Controlled responses variables 
f) = sideslip angle, deg 

p s — sin(atrim)p + cos(ortri m )r, roll rate about x stability axis, deg/s 
Control effectors 

Sdt — differential horizontal stabilizer, deg 
Sa = aileron, deg 
Sr— rudder, deg 

Srtv = differential pitch thrust vectoring, deg 
Syyv = yaw thrust vectoring, deg 
Actuator dynamics 
Sdt- 

30 2 a 

[0 707 30] rate ^ m * t = ^0 deg/s ampl. limit = ±17.5 deg 


75 2 

[0.59, 75] 


rate limit = 100 deg/s ampl. limit = ±27.5 deg 


rate limit = 100 deg/s ampl. limit = ±30.0 deg 


rate limit = 60 deg/s ampl. limit = ±30.0 deg 


Prefilters 


42613(1 ) 2 
(0)(100) 2 
22 . 8 ( 2 ) 

= ( 0 )( 100 ) 


G' P = G P - 


r> —r J 

G -- G "'W 


[0.707, 2] 


[0.707, 10] 


*{K (zi)/(pi)t?i, a>i3} = {tf (^ + zi )/(s + pi)[j 2 + + to 2 ]}. 

b Compensator prior to inclusion of s before limiter. 


Increasing the type q Lp of q Lp from 1 to 2 will improve the perfor- 
mance of the software limiters, provided, of course, that stability 
margins are not appreciably reduced by this increase, and they were 
not in this application. Note that the zeros in Eq. (29) were chosen 
to be a decade below the crossover frequencies of the associated 
loops so as to preserve the crossover characteristics of L p and L p . 
The final compensators are given in Table 1 along with the prefilter 
dynamics F p and F p . The bandwidths of the resulting closed-loop 
transfer functions, and pjp c , are determined by F p and F p . 
As in the case of the SISO example, in implementing the software 
limiters, the s appearing before each limiter was subsumed into the 
compensator that preceded it. 

Figures 10 and 1 1 show the p s and responses with and without 
the software limiters to inputs defined as follows. 

Figure 10: 

Pc(t) = o deg pc (t) = 180 • i(t) deg/s 



time, s 


Fig, 10 Roll rate response of system of Fig. 8a to input of Eqs. (30), with 
and without software rate limiter (rate and amplitude limited actuators 
always present). 



Fig. 11 Sideslip response of system of Fig. 8a to input of Eqs. (30) with 
and without software rate limiter (rate and amplitude limited actuators 
always present). 

Figure 11: 

&(0 = 30-i (0 deg p c {t) = 0 deg/s (30) 

where 

1(0 = «(0 — 2 u(t — 2) + 2 u(t — 4) — 2u(t — 6) + u(t — 8) 
and u{t) = unit step. 

As can be seen, the responses of the systems with software limiters 
are superior to those without. This is especially evident in the 
responses. Although not as dramatic a difference exists between 
the p s responses, the response of the system without the software 
limiters lags the response of the system with the limiters by as much 
as 0. 3-0.4 s for t > 4 s. In terms of handling qualities degradations, 
effective time delay increments of these magnitudes are certainly 
very serious. 15 It should be emphasized that in the simulation the 
actuators themselves always contain the rate and amplitude limiters 
shown in Fig. 9 and Table 1, and actuator amplitude limiting did 
occur with and without the software limiters. In addition, of course, 
software rate limiting occurred when these limiters were in place, 
and actuator rate limiting occurred when they were not. 

At the trim Mach number and altitude given in Table 1 , the large 
roll rates in Fig. 10 do not necessarily invalidate the use of a linear 
aerodynamic model. However, the large sideslip angles in Fig. 1 1 
for the system without software limiting are certainly in the non- 
linear aerodynamic range. Bearing in mind that this example was 
not intended to serve as an accurate predictor of vehicle response 
in large amplitude maneuvers but to demonstrate the performance 
improvements that can accrue with the proposed design approach, 
the results demonstrated in Figs. 10 and 1 1 were deemed acceptable. 
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V. Discussion 

The analyses of Secs. HI and IV have always begun with the 
assumption of closed-loop stability. Then necessary and sufficient 
conditions for asymptotic stability were used to derive a relation 
between the type of the injected signals D ni (s) in the quasilinear 
representation of the system and the type of the loop transmissions 
Lj ( 5 ) . Of course, stability under linear operation is guaranteed; how- 
ever, stability under limiting conditions is not. This is especially 
true if open-loop stability is not in evidence. However, when satu- 
ration occurs, and closed-loop stability is in evidence, the proposed 
technique will always result in improved performance as compared 
with a system with a similar feedback topology in which it is not 
employed. Stability limitations under saturation still need to be de- 
termined through simulation. 

Finally, it is imperative to point out that the rate limitations of 
flight control actuators are a function of aerodynamic loading. Thus, 
a practical implementation of the proposed design might require the 
elements of the K matrix in Eq. (14) and Fig. 8a to be scheduled 
with dynamic pressure. 

VI. Conclusions 

Based on the research just described, the following conclusions 
can be drawn. 

1) The inclusion of very simple software rate limiters in a flight 
control system with appropriate compensator and loop transmis- 
sion types may offer significant performance improvements in the 
presence of actuator rate saturation. 

2) As presented here, the design technique requires control ef- 
fectors to be ganged or grouped so that any effector or actuator is 
driven by a single compensated error signal. 

3) The design technique places restrictions upon the type of 
the loop transmissions and compensators, where type refers to the 
exponent on any free s in the denominator of the transfer func- 
tion. 

4) An approach was offered for determining the restrictions upon 
the transfer function type by examination of the steady-state behav- 
ior of the system. 

5) Although improved performance as compared with systems 
without the software limiters can be expected, no guarantees of 
closed-loop stability can be made. 
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